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ABSTRACT 

We investigate the relationship between the linewidths of broad Mg II A2800 and 
H/5 in active galactic nuclei (AGNs) to refine them as tools to estimate black hole (BH) 
masses. We perform a detailed spectral analysis of a large sample of AGNs at inter- 
mediate redshifts selected from the Sloan Digital Sky Survey, along with a smaller 
sample of archival ultraviolet spectra for nearby sources monitored with reverberation 
mapping (RM). Careful attention is devoted to accurate spectral decomposition, espe- 
cially in the treatment of narrow-line blending and Fe II contamination. We show that, 
contrary to popular belief, the velocity width of Mg II tends to be smaller than that 
of H/3, suggesting that the two species are not cospatial in the broad-line region. Us- 
ing these findings and recently updated BH mass measurements from RM, we present 
a new calibration of the empirical prescriptions for estimating virial BH masses for 
AGNs using the broad Mg II and H/3 lines. We show that the BH masses derived from 
our new formalisms show subtle but important differences compared to some of the 
mass estimators currently used in the literature. 
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1. Introduction 

It is generally accepted that active galactic nuclei (AGNs) are powered by the release of grav- 
itational energy from material accreted onto supermassive black holes (BHs). The determination 
of BH mass (Mbh) is crucial for understanding the AGN phenomena, the cosmological evolution 
of BHs, and even the coevolution of AGNs and their host galaxies. Yet, for such distant objects, it 
is currently impossible to obtain direct measurement of Mbh using spatially resolved stellar or gas 
kinematics. Fortunately, significant advances have been made in recent years from reverberation 
mapping (RM) studies of nearby Seyfert galaxies and quasi-stellar objects (QSOs; e.g., Wandel et 
al. 1999; Kaspi et al. 2000; Peterson et al. 2004). First, the anti-correlation between the radius 
of the broad-line region (BLR) and the velocity width of broad emission lines for single objects 
supports the idea that the BLR gas is virialized and that its velocity field is dominated by the grav- 
ity of the BH (Peterson & Wandel 1999, 2000; Onken & Peterson 2002). Second, the size of the 
BLR scales with the continuum luminosity (Kaspi et al. 2000, 2005), approximately as /? oc 
(Bentz et al. 2006, 2009); the R-L relation offers a highly efficient procedure for estimating the 
BLR size without carrying out time-consuming RM observations. And third, the BH masses esti- 
mated by RM are roughly consistent (Gebhardt et al. 2000b; Ferrarese et al. 2001; Nelson et al. 
2004; Onken et al. 2004) with the predictions from the tight correlation between Mbh and bulge 
stellar velocity dispersion established for inactive galaxies (the Mbh-i7* relation; Gebhardt et al. 
2000a; Ferrarese & Merritt 2000). These developments imply that we can estimate the BH mass in 
type 1 (broad-line, unobscured) AGNs by simple application of the virial theorem, Mbh = fRv^/ G, 
where / is a geometric factor of order unity that depends on the geometry and kinematics of the 
line-emitting region, R is the radius of the BLR derived from the AGN luminosity, and v is some 
measure of the virial velocity of the gas measured from single-epoch spectra. The feasibility of 
obtaining R and v from single-epoch spectra enables Mbh to be estimated very efficiently for large 
samples of AGNs, especially for luminous quasars at higher redshift that typically exhibit only 
slow and small- amplitude variability (e.g., Kaspi et al. 2007), with the assumption that the virial 
relation is independent of redshift and can be extrapolated to higher luminosities and masses. In 
practice, for those AGNs that have measurements of cr^, / is determined empirically by scaling the 
virial masses to the Mbh-c'"* relation of inactive galaxies (e.g., Onken et al. 2004). Implicit in this 
practice is the assumption — one open to debate (Greene & Ho 2006; Ho et al. 2008; Kim et al. 
2008) — that active and inactive BHs should follow the same Mbh-c* relation. The most widely 
used estimator for v is the full width at half-maximum (FWHM) of the line. 

Now, a large number of formalisms to estimate Mbh from single-epoch spectra have been pro- 
posed in the recent literature, using different broad emission lines optimized for different redshift 
regimes probed by (widely available) optical spectroscopy. At low redshifts, the lines of choice 
are H/3 (Kaspi et al. 2000; Collin et al. 2006; Vestergaard & Peterson 2006) or Ha (Greene & 
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Ho 2005). At intermediate redshifts, Mgll A2800 is used (McLure & Jarvis 2002), while at high 
redshifts, one has to resort to C IV A1549 (Vestergaard 2002; Vestergaard & Peterson 2006). These 
formalisms are ultimately calibrated against RM masses based on the H/? BLR radius and linewidth 
measured from the variable (rms) spectra (Peterson et al. 2004). Because the H/? linewidth is typi- 
cally smaller in the rms spectra than in the single-epoch or mean spectra (Vestergaard 2002; Collin 
et al. 2006; Sulentic et al. 2006), some authors have proposed that the FWHM used in the H/3- 
based formalisms should be further corrected to obtain unbiased Mbh estimates (Collin et al. 2006; 
Sulentic et al. 2006). 

As for the Mgll-based formalisms, because there are very few RM experiments of the Mgll 
line, they are either based on the RM data for H/9 (e.g., McLure & Jarvis 2002; McLure & Dunlop 
2004) or calibrated against the H/3 formalisms themselves (e.g., KoUmeier et al. 2006; Salviander 
et al. 2007). A strong, underlying assumption is that Mg 11 and H/3 are emitted from the same 
location in the BLR and have the same linewidth (see also Onken & KoUmeier 2008). In support 
of this assumption, some authors find that Mg II and H/? indeed have very similar linewidths (e.g., 
McLure & Jarvis 2002; McLure & Dunlop 2004; Shen et al. 2008; also cf. Salviander et al. 2007). 
However, there are conflicting results in the literature: Corbett et al. (2003) claimed that Mgll is 
generally broader than H/3, whereas Dietrich & Hamann (2004) came to an opposite conclusion. 
Certainly, the most direct way to settle this issue is through direct RM of the Mg II line. So far 
there are only two objects that have successful Mgn RM, NGC 5548 (Clavel et al. 1991; Dietrich 
& KoUatschny 1995) and NGC 4151 (Metzroth et al. 2006). These studies tentatively suggest that 
Mgll responds more slowly to continuum variations than H/3, implying that the MgU-emitting 
region is larger than that radiating H/3. 

Thus, there are still some important open questions regarding the robustness of Mbh measure- 
ments based on Mg II. What is the relation between the linewidths of H/3 and Mg II? Are estimates 
of Mbh based on Mg II consistent with those based on H/3? These basic questions are critical for 
understanding the systematic uncertainties in studies of the cosmological evolution of BHs (cf. 
Shen et al. 2008; McGill et al. 2008; Denney et al. 2009a). To address the above questions, we 
perform a detailed comparison of the widths of the Mgll and Hp lines using single-epoch spectra 
for a large, homogeneous sample of Seyfert 1 nuclei and QSOs at intermediate redshifts culled 
from the Sloan Digital Sky Survey (SDSS; York et al. 2000). We further compare single-epoch 
Mg II linewidths with H/3 linewidths measured from the rms spectra of AGNs with RM observa- 
tions, finding systematic deviations between the two. We present a recalibration of the Mg II virial 
mass estimator and compare our formalism with previous ones in the literature. 

This paper adopts the following set of cosmological parameters: Hq=70 km s"^ Mpc"^ , Q„=0.3, 
and nA=0.7. 
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2. Sample and Data Analysis 
2.1. The Samples 

The sample of best-studied H/3 emission lines is the one compiled by Peterson et al. (2004) 
for RM studies of 35 low-redshift AGNs. To compare H/3 and Mg II for this sample, we located 
usable ultraviolet (UV) spectra for 29 sources, 16 from the Hubble Space Telescope (HST) and 
13 from the International Ultraviolet Explorer (lUE) data archives. This sample will be used to 
study the relationship between single-epoch Mg II linewidths and H/? linewidths measured from 
rms spectra, and to fit a new Mbh formalism based on single-epoch Mg II. 

We also selected Seyfert 1 galaxies and QSOs in the redshift range 0.45 <z< 0.75 from the 
Fifth Data Release (DR5) of the SDSS spectroscopic database (Adelman-McCarthy et al. 2007). 
Within this redshift range, both H/? and Mgll lie within the SDSS spectral coverage. To ensure 
accurate measurement of both lines, we only select objects with a mean signal-to-noise ratio (S/N) 
> 20 per pixel in both the Hp (4600-5100 A) and the Mg II (2700-2900 A) regions. We discarded 
26 spectra that have either broad absorption lines or too many narrow absorption lines, or for 
which the H/3 or the Mg II regions were corrupted by bad pixels. The remaining 495 objects have 
spectra that can be well fitted, as confirmed by visual inspection. This sample will be used to 
investigate the FWHM relation between Mg II and H/3 in single-epoch spectra and to compare our 
Mg II formalism with others in the literature. 



2.2. Spectral Fitting 

The spectra are first corrected for Galactic extinction using the extinction map of Schlegel et 
al. (1998) and the reddening curve of Fitzpatrick (1999). Then the spectra are fitted using an IDL 
code based on MPFIT (Markwardt 2009), which performs -minimization by the Levenberg- 
Marquardt technique. Bad data are masked during the fitting. 

To measure the H/3 line, we perform continuum subtraction and emission-line fitting follow- 
ing the method described in detail in Dong et al. (2008). We first fit simultaneously the featureless 
nonstellar continuum (assumed to be a power law), the Fe II multiplet emission, and other emis- 

o 

sion lines in the wavelength range 4200-5600 A, giving emphasis on the proper determination of 
the local pseudocontinuum (continuum -i- Fe II emission). For spectra with fits having a reduced 
> 1-1 around H/? (4750-5050 A), a refined fit of the emission-line profiles is performed to the 
pseudocontinuum-subtracted spectra using the code described in Dong et al. (2005). Each line of 
the [Olll] AA4959,5007 doublet is modeled with two Gaussians, one accounting for the line core 
and the other for a possible blue wing as seen in many objects. The doublet lines are assumed to 
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have the same redshifts and profiles, and the flux ratio A5007/A4959 is fixed to the theoretical value 
of 3. The narrow component of H/3 is fitted with one Gaussian, assumed to have the same width 
as the line core of [O III] A5007. The broad component of H/? is fitted with as many Gaussians as 
statistically justified (see Dong et al. 2008 for details). 

To measure the Mg II line, we adopt the following procedure. We first obtain an initial estimate 
of the nonstellar featureless continuum by fitting a simple power law, 



to the data in several continuum windows near 2200, 3000, 4000, and 4200 A that suffer little from 
emission-line contamination, if available. Then, the power-law local continuum, a Balmer con- 
tinuum, and an Fell emission template, which together constitute the so-called pseudocontinuum, 
are fitted simultaneously. The fitting is performed in the restframe wavelength range 2200-3500 
A, if available, with the small region contaminated significantly by Mg II masked out. The fitting 
range is set by the wavelength coverage of the UV Fe II template, /^"^(A), which was generated 
by Tsuzuki et al. (2006) based on their measurements of IZw 1. In the wavelength region cov- 
ered by Mgll emission, they employed a semi-empirical iteration procedure to build the template. 
They first generated a theoretical Fe II model spectrum with the photoionization code CLOUDY 
(Ferland et al. 1998) and subtracted it from the observed IZw 1 spectrum around Mgll. Then the 
Mg II doublet was fit assuming each line has the same profile as Ha. And finally they obtained the 
Fe II template underneath Mg II by subtracting the Mg II fit from the observed spectrum. To match 
the linewidth and possible velocity shift of Fe II lines, we build the Fe II model by convolving the 
IZW 1 template with a Gaussian of width fig and shifting it with a velocity Vshift in logarithmic 
wavelength space (i.e., the velocity space because dlnA = dA/A = v/c), as follows. 



As in Dietrich et al. (2002), the Balmer continuum is assumed to be produced in partially 
optically thick clouds with a uniform temperature, Q 




(1) 



f\X;c,v,^ift,cT,) = c/i"^(A,v,hift)® G(A,ag) . 



(2) 



f'^iX;d,T„T,) = dBx{X,Te)il-e-^'); A < Abe 



(3) 




(4) 



We do not account for the velocity broadening of the Balmer continuum, because the Balmer continuum in our 
fitting range is insensitive to this effect. 
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where Abe = 3646 A (3.4 eV), Tbe is the optical depth at Abe, and 5a(A, Tg) is the Planck blackbody 
spectrum at the electron temperature T^. 

To sum up, the full model for the pseudocontinuum is as follows: 

/(A) = /L(A;a,/3)+/'=ii(A;c,VsMft,(Tg)+/«^^(A;J,r„r3H) . (5) 

The fitting is performed in logarithmic wavelength space. During the fitting, the normalization a 
and slope (3 of the power-law continuum, the normalization c, velocity shift VsMfi and broadening 
velocity tTg of the Fell emission, and the parameters d, Te, and Tbe of the Balmer continuum are set 
to be free parameters. 

We note that in the fitting range of 2200-3500 A the Balmer continuum is hard to be con- 
strained and separated from the power-law continuum and Fe II emission (cf . Figure 8 of Tsuzuki 
et al. 2006). In this work, we are not concerned with the properties of the Balmer continuum, but 
with the proper separation of the power-law continuum, Fe II, and Mg II. To minimize the effect 
of the possible poor fitting of Balmer continuum on the determination of the power-law contin- 
uum, Fe 11, and Mg II, we constrain the power-law continuum parameters in such a way that they 
vary only around the best-fit values obtained from the first step, by a factor of < 10% for the nor- 
malization a and < 20% for the slope /3. During this step, we assign additionally larger weights 
to the regions 2400-2650 A and 2920-2990 A in order to improve the fit for the Fe II emission 
surrounding Mgll (cf. Section 2.3 of Dong et al. 2008). 

Once the pseudocontinuum is fitted and subtracted, the Mg II emission line is fitted in the 
range of 2700-2900 A, if available. For the SDSS sample, there are a few cases where a small 
number of narrow absorption lines are present around the Mgll emission line. To further eliminate 
the absorption lines in these objects, we first fit the Mgll emission line with one Gaussian, and 
then mask those pixels of absorption features deviating strongly from the model. The Mg II line is 
fitted in the following way. Each of the two Mg II AA2796, 2803 doublet lines is modeled with two 
components, one broad and the other narrow. The broad component is a truncated five-parameter 
Gauss-Hermite series (van der Marel & Franx 1993; see also Salviander et al. 2007); the narrow 
component is a single Gaussian. The broad components of the doublet lines are set to have the 
same profile, with the flux ratio A2796/A2803 set to be between 2:1 and 1:1 (Laor et al. 1997), 
and the doublet separation set to the laboratory value. The same prescription is applied to the 
narrow components, with the following additional constraints: FWHM < 900 km s"^ and flux 
< 10% of the total Mgll flux (Wills et al. 1993; see also McLure & Dunlop 2004). The fitting 
results for all the 495 objects are reasonable according to our visual inspection finally. The FWHM 
value is measured from the Gauss-Hermite model of Mgll A2796. The monochromatic flux of the 
continuum is measured from the fitted power law. 

There are several other emission lines in the fitting region, identified from the composite 
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SDSS QSO spectrum (see Table 2 of Vanden Berk et al. 2001); yet, because of their weakness, we 
simply masked them out in the fit. Because of the limited wavelength coverage of the RM sample, 
we cannot separate the Balmer continuum from the power-law continuum. Thus, the Balmer con- 
tinuum was not included in the fits for this sample. Additionally, there are deep narrow absorption 
features around Mgll in the spectra of NGC 3227, NGC 3516, NGC 3783, and NGC 4151 in the 
RM sample. Each of the absorption features is fitted simultaneously with a Gaussian when fitting 
the Mg II emission line. 

We estimate the measurement uncertainties of the parameters using the bootstrap methocj^ 
described in Dong et al. (2008, Section 2.5). The estimated 1 a errors for the broad-line fluxes 
are typically 10% for Mgll and 8% for H/?, while the errors on the broad-line FWHM are ~20% 
for Mgll and ~15% for The power-law continua have uncertainties of 8% for the slope and 
5% for the normalization. The above discussion does not account for possible systematic errors 
resulting from the subtraction of the continuum or our treatment of the Fe II and narrow lines. 

Figure 1 shows two examples of the fits. The continuum and emission-line parameters for the 
RM and SDSS samples are listed in Tables 1 and 2, respectively. The data and fitting parameters 
are available online for the decomposed spectral components (continuum, Fe II, and other emission 
lines). 



2.3. Regression Methods 

In the next section, we will fit the linewidth relations and BH mass estimators using several 
regression methods. The purpose of using these different methods is for ease of comparison with 
results in the literature, and to investigate possible differences in the fitting results caused by the 
different methods. Here we briefly summarize the regression methods used. 



-To estimate the errors on the fitted parameters, we generate 500 spectra by randomly combining the scaled model 
emission lines of one object (denoted as "A") to the emission-Une subtracted spectrum of another object (denoted as 
"B"). The emission-line model of object "A" is scaled in such a way that it has the same broad Mg II flux as object "B," 
in order to minimize changes in S/N within the emission-line spectral regions in the simulated spectra. Then, we fit 
the simulated spectra following the same procedure as described in Section 2.2. For each parameter, we consider the 
error typical of our sample to be the standard deviation of the relative difference between the input and the recovered 
parameter values. These relative differences turn out to be normally distributed for each of the parameters concerned. 

^Available at 

|http://staff.ustc. edu.cn/~xbdong/Data_Release/MgII_Hbeta/ together with auxiliary code to explain the parameters 
and to demonstrate the fitting. 
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1. Ordinary least-squares (OLS), which is a least-squares regression method without consider- 
ing measurement errors. 

2. Weighted least-squares (WLS), which takes into account only measurement uncertainties in 
the dependent variable. 

3. FITexy (Press et al. 1992), which numerically solves for the minimum orthogonal us- 
ing an interactive root-finding algorithm. It accounts for measurement uncertainties in both 
coordinates, but does not account for intrinsic scatter. 

4. FITexy _T02, the version of FITexy modified by Tremaine et al. (2002), accounts for possible 
intrinsic scatter in the dependent variable by adding in quadrature a constant to the error value 
so as to obtain a reduced of 1. 

5. Gaussfit (McArthur et al. 1994), which implements generalized least-squares using the ro- 
bust Householder Orthogonal Transformations (Jefferys 1980, 1981). It can handle errors in 
both coordinates, but does not account for intrinsic scatter. 

6. The bivariate correlated errors and intrinsic scatter (BCES) regression method (Akritas & 
Bershady 1996), which accounts for measurement errors on both coordinates in the fit using 
bivariate correlated errors, and possible intrinsic scatter (but does not output any quantifica- 
tion of this scatter). The results of the two symmetrical versions, bisector and orthogonal, 
are used in this paper. 

7. LINMIX_ERR (Kelly 2007), which accounts for measurement errors, nondetections, and 
intrinsic scatter by adopting a Bayesian approach to compute the posterior probability dis- 
tribution of parameters, given observed data. We also consider the multivariate extension, 
MLINMIX_ERR. 

As we find below, most of the above regression methods give consistent results. For the 
linewidth-linewidth relationships, since there is no prior knowledge about which variable is in- 
dependent and which is dependent, we adopt the results given by the BCES orthogonal method, 
which treats both variables symmetrically. For the BH mass scaling relations, we adopt formally 
the results given by the LINMIX_ERR method, since it is argued to be among the most robust 
regression methods with the possibility of reliable estimation of intrinsic dispersion (Kelly 2007). 

For some of the regression methods listed above no intrinsic scatter ((Jint) can be inferred. 
We can give a rough yet simple estimate of a-yai by deducting the contribution of the measure- 
ment errors from the variance in the regression residuals (cTtot^), by using an approximate rela- 
tion (T?t = cTtot" (cm)^, where ((Tm) is the median of the total measurement errors computed from 
Cm = (T3,-i-slope*(T;c. If (Ttot is Smaller than ((Tm), (Tmt was set to be 0. As a check, this rough estimate 
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can be compared with the intrinsic scatter given by some of the regression methods that provide 
such a measure. 



3. Results 



The main motivation of this work is to investigate whether reliable BH masses can be esti- 
mated using the Mgll linewidth as a virial velocity indicator. Linewidths are commonly parame- 
terized as FWHM, or sometimes as — the line dispersion or second moment of the line profile 
(Peterson et al. 2004). Both quantities have intrinsic strengths and weaknesses (see Section 3 of 
Peterson et al. 2004). Collin et al. (2006), in particular, argued that the use of FWHM rather 
than cTiine introduces systematic bias in Mbh estimates. Although aiine is a better tracer of virial 
velocity than FWHM in rms spectra (Peterson et al. 2004), the line dispersion is very sensitive to 
measurement errors in the line wings, making it especially susceptible to inaccuracies caused by 
deblending and subtraction of Fell and other emission lines, effects that are particularly significant 
in mean and single-epoch spectra. By contrast, the FWHM is less prone to these effects; it is more 
sensitive to corrections for the narrow-line component, which, fortunately, is quite weak for Mg II 
(see Section 4.1). In this work, we opt to use the FWHM to parameterize the linewidth. 



3.1. Single-epoch Mg n FWHM versus H/3 FWHM 

We first investigate the relation between the FWHM of Mg II and H/3, using single-epoch data 
from our SDSS sample. The relation is illustrated in Figure 2. A strong correlation is present, 
but apparently deviates from one-to-one. This trend has been noticed in the literature, but it was 
less prominent because of the narrower dynamical range in velocity covered in previous studies 
(e.g., Salviander et al. 2007; Hu et al. 2008). With our high-quality data, we can now fit a strict 
relation. We perform a linear regression in log-log space using the methods described in Section 
2.3; the results are listed in Table 3. As can be seen, most of the methods give mutually consistent 
results. For our subsequent analysis, we adopt the BCES (orthogonal) method because it treats 
both variables symmetrically (Section 2.3). We find 



log 



FWHM(Mgll) 
1000 km s-i 



= (0.81 ±0.02) log 
+ (0.05 ±0.01). 



FWHM(H/3) 
1000 km s-i 



(6) 



This means that the line-emitting locations of HP and Mgll in the BLR are not identical. If they 
were, we would expect a linear relation between the two, with no offset. The intrinsic scatter of 
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this relation, as given by the regression methods listed above, is extremely small and negligible 
compared to the measurement errors. The latter is actually comparable to the total scatter {a^ox) of 
the relationship, which is found to be 0.08 dex. 



3.2. Single-epoch Mg II FWHM versus rms H/3 a]^^ 

Since the assumption that Mg 11 FWHM is identical to H/3 FWHM does not hold, we explore 
the relation between Mg II FWHM and rms H/3 crime, which has been argued to be a good tracer of 
the virial velocity of the BLR clouds emitting (variable) H/3 (see references in Section 1). We use 

data for the 29 objects in the RM sample that have UV spectra to perform this exploration. Mg II 
FWHM is measured from the single-epoch HST/IUE spectra, as listed in Table 1 . The data for 
rms H/? o"iine are mainly taken from Peterson et al. (2004). In addition, we use updated RM data 
for NGC4051 (Denney et al. 2009b), NGC4151 (Metzroth et al. 2006), NGC4593 (Denney et 
al. 2006), NGC5548 (Bentz et al. 2007), and PG 2130+099 (Grier et al. 2008). For objects with 
multiple measurements, the geometric mean (i.e., the mean in the log scale) was used. 

We find that the slope of the relation between Mg 11 FWHM and rms H/3 crime deviates from 
unity, with a best fit of 



log 



q-iine(H/j,rms) 
1000 km s-i 



= (0.85±0.21)log 
- (0.21 ±0.12) . 



FWHM(Mgll) 
1000 km s-i 



(7) 



The formal relation is nonlinear although the significance level is only about 1 cr. A nonlinear 
relation between Mg n FWHM and rms H/3 a\{-as is not very surprising, in light of a similar situation 
observed for H/3 FWHM (Collin et al. 2006; also Sulentic et al. 2006, Section 1). For verification, 
we also fit the relation between H/9 FWHM in the mean spectra and rms H/? aune using data for 
35 objects in the RM sample; the FWHM data are taken from Collin et al. (2006) and from the 
updated sources mentioned above. The best-fit relation deviates from unity even more seriously 
than the case of Mg II FWHM: 



log 



criine(H/j,rms) 
1000 km s-i 



= (0.54 ±0.08) log 
- (0.09 ±0.05) . 



FWHM(H/3,mean) 
1000 km s-i 



(8) 



These relations between rms H/3 crime and Mgn and H/? FWHM are illustrated in Figure 3. The 
total 1 cr scatter around these relationships is cr,ot = 0.12 dex for Equation (7) and 0.09 dex for 
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Equation (8). Given the relatively small measurement errors of the linewidths, there likely exists 
intrinsic scatter in these relationships. Using the simple method of deducting the measurement 
errors from the total scatter, as described in Section 2.3, we find a"int~ 0.09 dex and ~ 0.08 dex for 
the relationships of Equations (7) and (8), respectively. The underlying reason for the nonlinearity 
of these relationships may be, at least partially, that the aiine of rms spectra traces the velocity of 
the line-emitting region that responds to continuum variation, while the FWHM of single-epoch 
spectra may be contributed by various components (see Section 4.2 for a discussion). 



As described above, Mg II FWHM is not identical to, but rather generally smaller than, H/3 
FWHM; for Mgll FWHM > 6000 km s"\ the difference is > 0.2 dex. This means that one of the 
fundamental premises of the previous Mg ll-based formalisms — that Mg II and H/? trace similar 
kinematics — does not hold. Moreover, similar to the behavior of H/3 FWHM, Mg II FWHM seems 
not to be linearly proportional to rms H/3 crune • If nns crune is more directly linked to the virial 
velocity, this implies that we cannot build a virial Mbh formalism by simply assuming Mbh oc 
FWHM^. Furthermore, the Mbh data of the RM AGNs used in McLure & Jarvis (2002) and 
McLure & Dunlop (2004) have since been recalibrated or updated (Peterson et al. 2004; Denney 
et al. 2006, 2009b; Metzroth et al. 2006; Bentz et al. 2007; Grier et al. 2008). Thus, it is necessary 
to reformulate the virial Mbh formalism based on single-epoch Mgll FWHM. 

We proceed by assuming that there is a tight relation between the BLR radius of the Mg II- 
emitting region and the AGN continuum luminosity, in the form Rugii oc L^, and another between 
the virial velocity of Mgll and the FWHM of the line, in the form v^j^j^^i oc FWHM^. Then, using 
the 29 objects with the Mbh values based on RM and the Mgll data measured here (Table 1), we 
calculate the free parameters by fitting 



where L3000 = AL;v(3000 A). The RM-based Mbh data are mainly taken from Peterson et al. (2004), 
who calibrated the /-factor by normalizing to the Mbh-c"* relation of Onken et al. (2004); Mbh for 
the updated objects comes from the references given in Section 3.2. 



3.3. Practical Formalism for New Mg n-based Mbh Estimator 




+ 7 log 



FWHM(Mgll) 



(9) 



1000 km s-i 



We fit Equation (9) following four schemes, using the (LINMIX_ERR/MLINMIX_ERR) 
method of Kelly (2007): 
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1. a, (3, and 7 are treated as free parameters. 

2. a and /3 are treated as free parameters, but, as in all previous formalisms, we fix 7 = 2. 

3. a and /? are treated as free parameters, but we set 7 = 1.70, as suggested by Equation (7). 

4. a and 7 are treated as free parameters, but we fix /3 = 0.5, as suggested by the latest R—L 
relation (Bentz et al. 2006, 2009). 

Table 4 lists the best-fit regression for each scheme (Columns 1^), as well as comparisons 
between the Mbh estimates based on each scheme and the RM-based masses (Column 5). It is 
apparent that the best-fit values for P for all the schemes are consistent with 0.5 within 1 a error. 
Interestingly, 7 appears to be marginally smaller than 2, since the standard deviation of the BH 
mass for Scheme 2 is slightly larger than that for the other three schemes. If we set P - 0.5 (i.e., 
adopt Scheme 4), the best- fit Mg ll-based formalism is 



log 



/ M 



BH 



= (1.13±0.27) + 0.51og(^^ 
+ (1.51 ±0.49) log 



'3000 



erg s" 
FWHM(Mgll) 



1000 km s- 



(10) 



Fitting the H/3 FWHM data for the 35 RM objects under the same assumptions (/3 = 0.5), the 
H/3-based formalism becomes 



log 



BH 



V 106 M, 







(1.39±0.14)+0.51ogf ^ 

\10^^ erg s ^ J 



+ (1.09 ±0.23) log 



FWHM(H/5) 
1000 km s-i 



(11) 



The best-fitting 7 = 1.09 ±0.23 agrees well with the (7ii„e(H/3,rms)-FWHM(H/3) relation derived 
in Equation (8). 

Comparisons between the RM-based masses and the Mbh estimates from our new Mbh for- 
maUsms using Mgll FWHM (Equation (10)) and H/3 FWHM (Equation (11)) are illustrated in 
Figure 4. Following Vestergaard & Peterson (2006), we calculate the deviation of the new cal- 
ibrated single-epoch Mbh estimates from the RM-based masses, AMbh(RM). The mean of the 
deviations, (AMbh(RM)), is only 0.01 dex for our Mgll estimator, and the 1 a scatter is 0.4 dex 
(Column 5 in Table 4). As a comparison, if we use the formalism of McLure & Dunlop (2004), 
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the deviations from the same RM-based masses have a mean of 0.38 dex and a 1 a scatter of 0.45 
dex. For our H/3 estimator, (AMbh(RM)) is 0.01 dex and the 1 a scatter is 0.3 dex, compared 
to ( AMrh(RM)) = 0.05 dex and a = 0.4 dex if the formalism of Vestergaard & Peterson (2006) 
is usedo It should be noted that these scatters of the scaling relations, which give a measure of 
the uncertainty in estimating Mbh from the single-epoch spectroscopic data, is relative to the RM- 
based masses only. Since, as pointed out by Vestergaard & Peterson (2006), the RM-based masses 
themselves are uncertain typically by a factor of ~ 2.9 (as calibrated against the Mbh-o"* relation; 
Onken et al. 2004), the absolute uncertainty of the masses thus estimated is even higher. For the 
H/3 formalism, we find this absolute uncertainty to be a factor of ~ 3.5, to be compared with a 
factor of ~ 4 given in Vestergaard & Peterson (2006); for Mgll, we estimate that the absolute 
uncertainty is a factor of ~ 4. 

As shown above, our new formalisms improve somewhat the scatter in the single-epoch Mbh 
estimates compared to previous H/? and Mgll estimators, by 0.1 dex and 0.05 dex, respectively. 
Given the same linewidth and luminosity data used in this work and in Vestergaard & Peterson 
(2006), the reduction in the scatter of the A/bh estimates should result from a decrease in the intrin- 
sic dispersion of our improved single-epoch Mbh formalisms. Using the LINMIX_ERR method, 
the intrinsic scatter inherent in our Mbh formalism can be inferred to be 0.08 dex (1 cr) for the H/3 
and 0.14 dex for Mg II. 

Figure 5 compares our new Mg Il-based formalism (we show only Schemes 2 and 4) with the 
previous H/5-based formalisms of Vestergaard & Peterson (2006; panel a) and Collin et al. (2006; 
panel b), as well as our newly derived version using the SDSS sample of 495 Seyfert Is and QSOs 
(Equation (11); panel c). The Mbh residuals between our Mg II formalism and the H/? formalisms 
are listed in Table 4 (Columns 6-8). While our Mgll-based formalism, especially for Scheme 4 
(Equation (10)), agrees well with our H/3-based formalism (Equation (11)), note that it deviates 
markedly from the H/3 formalism of Vestergaard & Peterson. This confirms previous suspicions 
(Collin et al. 2006; Sulentic et al. 2006) that the use of H/3 FWHM from mean and single-epoch 
spectra with the assumption 7 = 2 introduces systematic bias into Mbh estimates. 

We further compare our new Mgll-based formalism (Equation (10)) with other Mgll for- 
malisms widely used in the literature. Figure 5 illustrates that the Mbh estimates following the 
formalisms of McLure & Dunlop (2004; panel d), KoUmeier et al. (2006; panel e), and Salviander 
et al. (2007; panel f) show large systematic deviations, mostly in the sense of being smaller than 
ours. The deviations stem primarily from the recalibration of the RM masses; other factors are 



It should be noted that here we use the averaged spectrum for an object with more than one observation, unlike 
in Vestergaard & Peterson (2006) where the individual single-epoch H/3 spectral data were used in the regression. 
If we take the latter approach, our H/3 estimator gives a (AMbh(RM)) = -0.07 dex and a scatter of 0.33 dex, while 
Vestergaard & Peterson (2006) gave (AMbh(RM)) = -0. 12 dex and a scatter of 0.45 dex. 
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discussed in Section 4.3. We note that a yet-unpublished Mg ll-based formalism by M. Vestergaard 
et al. (in preparation) used in the recent literature (e.g., Kelly et al. 2009) is almost identical to our 
Scheme 2 (with 7 fixed to 2). 



4. Discussion 

4.1. Testing the Effect of Narrow-line Subtraction 

The narrow component of Mgll is generally weak in luminous type 1 AGNs (e.g., Wills et 
al. 1993; Laor et al. 1994), and so its contribution to the total line flux can be safely neglected. 
However, its presence might have a more pronounced impact on the FWHM measurement of broad 
Mgll. In the literature, narrow Mgll was accounted for in the line fitting by some authors (e.g., 
McLure & Dunlop 2004), but not by others (e.g., Salviander et al. 2007). As there is usually no 
clear inflection in the Mg II profile, separating narrow Mg U from the broad component is often 
challenging. Fortunately, for the spectra in our SDSS sample, [O III] A5007 is present, and thus 
we can use [O III] to try to constrain narrow Mg 11, to test the effect of narrow Mg 11 on the FWHM 
measurement of broad Mg II, and also to test the reliability of our Mg II fitting strategy. 

In addition to the default fitting strategy described in Section 2.2, in which narrow Mgll 
is modeled as a single free Gaussian, we tried two alternative strategies in which narrow Mg II 
is (A) not fit at all, and (B) is fit using a single-Gaussian model constrained to that of the line 
core of [OlII]. Broad Mgll is modeled as described in Section 2.2. We find that, for the 495 
objects in our SDSS sample, the distributions of the reduced of the Mg 11 emission line fit of 
the three approaches can be approximated reasonably well with a log-normal function. The peak 
and standard deviation of the reduced are very similar for all three, being (0.97, 0.10 dex) for 
the default strategy, (1.01, 0.10 dex) for Strategy A, and (0.99, 0.10 dex) for Strategy B. Regarding 
the FWHM of broad Mgll, the mean and standard deviation are (-0.04, 0.05) for log[ p^™4a^^^ ] 
and (0.00, 0.05) for log[ p^™4!uit) ]- Strategy B, the fitted flux of narrow Mgn is less than 
10% of the total line flux for almost all the objects. Our tests show that omitting the subtraction of 
narrow Mg n has a negligible effect on the FWHM of broad Mg II, typically decreasing it only by 
a tiny factor of 0.04 dex. We further confirm that our default procedure for modeling narrow Mg II 
is consistent with that using the [O III] core as a template. 

The above Strategy A is exactly the same as the Mg ll-fitting method adopted by Salviander et 
al. (2007). We also compared our method with that of McLure & Dunlop (2004). When fitting the 
spectra in our SDSS sample by the method of McLure & Dunlop (2004), on average the FWHM 
of broad Mg II is larger than that of our method by 0. 1 dex. 
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4.2. Mbh Estimators with Single-epoch H/3 and Mgn 

As analyzed in detail by Sulentic et al. (2006), the overall profile of the H/3 emission line, 
as viewed in single-epoch spectra, likely comprises multiple components emitted from different 
sites. First, as a recombination line, H/3 can arise from BLR gas that is very close to the central 
engine. Then H/3 can be gravitationally redshifted, as (part of) the component of the "very BLR" 
(Marziani & Sulentic 1993). Such clouds may be optically thin to the ionizing continuum, such 
that H/3 is no longer responsive to continuum variation (Shields et al. 1995). Second, like CiV 
A 1549, H/3 can be produced partly in high-ionization winds, as some observations suggest (see 
Marziani et al. 2008, and references therein). This wind component would not be virial. Third, H/3 
can also be produced on the surface of the accretion disk, both by recombination and coUisional 
excitation (Chen & Halpem 1989; Wang et al. 2005; Wu et al. 2008); this component would be 
highly anisotropic (cf. Collin et al. 2006). Considering the above factors, it is not surprising that 
single-epoch FWHM is not linearly proportional to crime for rms H/3 (Equation (8)). 

Mg II, as a low-ionization, coUisionally excited emission line, cannot be produced in clouds 
very close to the central engine. Furthermore, because Mg II originates only from optically thick 
clouds, radiation pressure force cannot act on them very significantly (cf. Marconi et al. 2008, 
2009; Dong et al. 2009a,b), and thus Mgll suffers little from nonvirial motion. Hence, compared 
to H/3, the FWHM of single-epoch Mg II should, in principle, deviate less, if at all, from the true 
virial velocity of the line-emitting clouds. This is suggested by the best-fit value for 7 in Equation 
(10), which indicates v^j^^ oc FWHM(MgII)i-^i^°-^^ 

Previously, researchers have feared that the substantial contamination of the Mg II region by 
Fe II multiplets might introduce significant uncertainties in its linewidth measurements, such that 
Mg ll-based Mbh estimates may not be as accurate as those based on H/3. With the recent avail- 
ability of a more refined UV Fe II template (Tsuzuki et al. 2006; cf. Vestergaard & Wilkes 2001), 
we have higher confidence that the linewidth measurements of Mg II are reasonably robust. Our 
work suggests that we can measure Mgll FWHM typically to within an uncertainty of ~20%. 
Nevertheless, it would be highly desirable to attempt to further improve the methodology for Fe II 
subtraction, not only in the UV but also at optical wavelengths. 



4.3. Comparison with Previous Studies 

As shown in Sections (see Figure 5), our Mgll- and H/3-based Mbh formalisms show, in 
addition to somewhat improved internal scatter, subtle but systematic deviations from some of 
the commonly used Mbh estimators in the literature. In general, the formalism prescribed by 
our Scheme 4 (Equation (10); Mbh oc FWHM^'^^^^ "^^) gives progressively higher and lower Mbh 
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values toward the low- and high-mass ends, respectively. The only exception is the H/3-based 
formalism of Collin et al. (2006), which gives roughly consistent results as ours over a relatively 
large mass range. The discrepancies between previous mass estimators and ours arise from one, 
or a combination, of the following factors incorporated into our analysis. (1) We use the most 
recently recalibrated and updated RM Mbh measurements from the literature (Peterson et al. 2004; 
Denney et al. 2006, 2009b; Metzroth et al. 2006; Bentz et al. 2007; Grier et al. 2008). (2) Our 
new formalism (Scheme 4, Equation (10)) uses the best-fitting value of 7 instead of the canonical 
value of 7 = 2. (3) For Mg II, we determine the scaling factor (incorporated into the coefficient 
a of Equation (9)) and the power-law index (/?) of the i?Mgii-L relation by fitting Equation (9) to 
the data, instead of simply using the existing /?h/3-L relation as a surrogate. (4) Differences in the 
line-fitting and determination of the FWHM. We discuss each of these factors in detail below. 

Specifically, assuming the canonical value of 7 = 2 in Equation (9) would, compared to our 
Scheme 4, underestimate Mbh at the low-end and overestimate Mbh at the high-end, for both 
Mg II and H/3. This accounts for most of the deviations from Vestergaard & Peterson (2006) and 
KoUmeier et al. (2006), and partially from others in Figure 5. In order to account for systematic 
biases with respect to RM-based Mbh, Collin et al. (2006) introduced a correction factor, which is 
dependent on H/9 FWHM, into their H/3-based formalism assuming 7 = 2. This correction has a 
similar effect as fitting 7 as a free parameter, as we do here, and thus the rough consistency between 
our results and theirs is not surprising. Factor (3) is also partially responsible for producing the 
deviations from some of the previous Mg Il-based Mbh estimators, such as those of McLure & 
Dunlop (2004) and Salviander et al. (2007), who assumed that both Mg n and H/3 obey the same 
R-L relation, and of KoUmeier et al. (2006), who used a very steep relation of i?Mgii L°^^. 

There have been previous reports of discrepancies between Mg n- and H/?-based estimators, 
which are sometimes claimed to correlate with luminosity or Eddington ratio (e.g., KoUmeier et al. 
2006; Onken & KoUmeier 2008). These effects can be traced, at least partially, to the one-to-one 
relation assumed between FWHM(Mg n) and FWHM(H/3), which is contradictory to the nonlinear 

relation found in this work. In fact, by adopting a nonlinear FWHM(Mg II)-FWHM(H/3) relation 
and Rugii oc LP'^, our new Mgll- and H/5-based estimators yield mutually consistent results for 
the SDSS sample (Figure 5, panel c). We verified that the previously claimed correlations of the 
residuals of the Mg II- and H/3-based estimators with luminosity or Eddington ratio largely vanish; 
a Spearman rank analysis indicates a chance probability of 0.09 for the former and 0.05 for the 
latter. 

It is generally accepted that the width of the variable part of the line, the line dispersion 
(Jiine measured in the rms spectrum, is by far the best tracer of the virial velocity of the BLR gas 
responsible for the variable portion of the emission line (e.g., Onken & Peterson 2002), such that 
Mbh (X c^jjg according to the virial theorem. If the virial velocity is estimated using FWHM (or 
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any other measure of the Unewidth) in single-epoch spectra, as long as its relation with a]ine is 
nonlinear, 7 in Equation (9) is expected to deviate from 7 = 2. This is exactly what we find in this 
work for FWHM(Mg II) at 1 a significance level, as well as for FWHM(H/5) at 6 a significance 
level. In fact, the fitted value of 7 = 1.51 ±0.49 for Mgll (Equation (10)) and 1.09 ±0.23 for H/9 
(Equation (11)) are almost identical to those derived from the fitted o'iine-FWHM(Mg ll) and aum- 
FWHM(H/3) relations, which have slopes of 1.70 ±0.42 and 1.08 ±0.16, respectively. Factor (3) 
is justified by the compelling evidence presented in this work that the line-emitting locations of 
up and Mg II in the BLR are not identical. Possible physical processes underlying factors (2) and 
(3) are discussed in Section 4.2. 

As an additional consideration, we have performed in this work refined and careful line fitting 
and determination of the FWHM, which may have subtle differences from previous results. These 
differences may also give rise to, to some extent, the systematic discrepancies between the mass 
relations in Figure 5 since we use our measured FWHM and luminosity data when producing the 
figure. For example, Salviander et al. (2007) did not subtract the narrow component of Mgll, 
leading to Mgll FWHM statistically smaller than ours; so the true deviations in Figure 5 (panel f) 
would be larger if their FWHM data were used. On the contrary, McLure & Dunlop (2004) over- 
subtracted the narrow component compared to ours (they considered a possible narrow component 
as having an upper limit of FWHM = 2000 km s"^ , much larger than the 900 km s"^ used in our 
work), and their Mgn FWHM are statistically larger than ours; thus, the true deviations in Figure 
5 (panel e) would be somewhat smaller if their FWHM data were used. 

Finally, the appropriateness of our approach is further justified by the fact that our formalisms 
give Mbh values consistent with the RM measurements with the least systematic bias, as well as 
a reduced (intrinsic) scatter compared to previous formalisms (Section 3.3). Moreover, we find 
consistent masses between the Mgll- and H/3-based estimators. We thus conclude that our Mbh 
estimators introduce less systematic bias compared to previous formalisms. Obviously, more RM 
measurements (for both H/3 and Mg n) are needed in order to improve the determination of the 
(7iine-FWHM relation, the R-L relation, and the index 7 in the Mbh FWHM''' relation. 



5. Summary and Conclusions 

We investigate the relation between the velocity widths for the broad Mg II and H/3 emission 
lines, derived from FWHM measurements of single-epoch spectra from a homogeneous sample of 
495 SDSS Seyfert Is and QSOs at 0.45 <z< 0.75. Careful attention is devoted to accurate spec- 
tral decomposition, especially in the treatment of narrow-line blending and Fe II contamination. We 
find that Mg II FWHM is systematic smaller than H/? FWHM, such that FWHM(Mg II) oc FWHM(H/5)' 
Using 29 AGNs that have optical RM data and usable archival UV spectra, we then investigate 
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the relation between single-epoch Mg II FWHM and rms 11(3 (line dispersion), a quantity 

regarded as a good tracer of the virial velocity of the BLR clouds emitting the variable H/9 com- 
ponent. We find that, similar to the situation for the FWHM of single-epoch H/3, single-epoch 
Mgll FWHM is unlikely to be linearly proportional to rms H/5 anne ■ The above two findings 
suggest that a major assumption of previous Mgll-based virial BH mass formalisms — that the 
Mg Il-emitting region is identical to that of H(3 — is problematic. This finding and the recent up- 
dates of the reverberation-mapped BH masses (Peterson et al. 2004; Denney et al. 2006, 2009b; 
Metzroth et al. 2006; Bentz et al. 2007; Grier et al. 2008) motivated us to recalibrate the Mbh 
estimator based on single-epoch Mg II spectra. 

Starting with the empirically well-motivated BLR radius-luminosity relation and the virial 
theorem, Mbh oc L'^FWHM''^, we fit the reverberation-mapped objects in a variety of different ways 
to constrain (3 and 7. For all the strategies we have considered, (3 has a well-defined value of ~0.5, 
in excellent agreement with the latest BLR radius-luminosity relation (Bentz et al. 2006, 2009), 
whereas 7 ~ 1.5 ±0.5, which is marginally in conflict with the canonical value of 7 = 2 normally as- 
sumed in past studies. Performing a similar exercise for H/? yields Mbh oc L°^FWHM(H/3)' °^^°'^^, 
which again significantly departs from the functional forms used in the literature. The 1 a uncer- 
tainty (scatter) is of the order of 0.3 dex relative to the RM-based masses for the H/3 estimator, 
and ~0.4 dex for the Mgll estimator. Using the same data set, the scatter of our H/3 mass scaling 
relation is reduced by 0.1 dex over that of Vestergaard & Peterson (2006), indicating improvement 
in the internal scatter. 

We use the SDSS database to compare our new Mbh estimators with various existing for- 
malisms based on single-epoch H/5 and Mg II spectra. BH masses derived from our Mg Il-based 
mass estimator show subtle but important deviations from many of the commonly used Mbh esti- 
mators in the literature. Most of the differences stem from the recent recalibration of the masses 
derived from RM. Researchers should exercise caution in selecting the most up-to-date Mbh esti- 
mators, which are presented here. 
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Fig. 1. — Examples of Mgll fitting for (top) the HST spectrum of Mrk 335 and (bottom) the 
SDSS spectrum of SDSS J13 1459.75+505932.9. The data are shown in black, the power-law 
AGN continuum in blue, the pseudocontinuum (power law plus Fell emission) in pink, the final 
model for all fitted components in green, and the continuum- subtracted emission-line spectrum in 
gray. For the multi-Gaussian fit to Mg II, the narrow components are shown in navy, the individual 
broad components in brown, the sum of all the broad components in cyan, and the total model 
(narrow plus broad) in red. 
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Fig. 2.— FWHM(Mgll) vs. FWHM(H/3) for our SDSS sample. The solid line represents the 
best- fitting power law with index 0.81. The dashed line represents a 1:1 relationship. A typical 1 
a error bar is also shown (top-left). 
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Fig. 3. — o-iine(H/3,rms) vs. FWHM(Mg II) of the 29 objects with Mgll FWHM measured in the 
paper (left panel) and FWHM(H/3, mean) of the 35 objects from Collin et al. (2006) and recent 
updated data (right panel). The solid lines show the best-fitting relations. The error bars are at 1 a. 
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Fig. 4. — BH masses estimated from RM plotted against masses obtained from Mgll (left panel; 
using the 29 objects with Mg II data measured in the paper), and from H/5 (right panel; using the 
35 objects from Peterson et al. (2004) and recent updated data). The solid line represents a 1:1 
relationship. The error bars are at 1 a. 
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Fig. 5. — Comparison of Mbh estimated by our Mgll formalism with other formulae discussed 
in Section 3, using the SDSS sample. The ^-coordinates of the green and red points represent 
Mgll masses estimated by Scheme 2 (Mbh oc FWHM^) and Scheme 4 (Mbh oc FWHM^-^^=^°'*''), 
respectively. The inset in each panel plots the histograms of Mbh; green and red lines denote the 
Mgll-based masses from our Schemes 2 and 4, and black lines are the comparison masses. The 
top three panels compare our Mgll masses with masses derived from different H/3 formalisms. 
The H/3 masses from Vestergaard & Peterson (2006; a) are systematically different from our Mg II 
masses, while those of Collin et al. (2006; b) are roughly consistent, and the best agreement 
comes from our newly derived formalism (Equation (9); c, red points). The bottom three panels 
compare our Mg Il-based masses with previous Mg Il-based formalisms: McLure & Dunlop (2004; 
d), KoUmeier et al. (2006; e), and Salviander et al. (2007; f). All show systematic deviations, 
mostly in the sense of giving lower masses than our formalism. 
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Table 1. Data for the Reverberation-mapped Sample 



Name 


DatalD 


logLsooo 


FWHM(Mgn) 


logLsioo 


FWHM(H/3, mean) 


o-;,„e(H/3, rms) 


logMBH(RM) 






(erg s"') 


(kms-i) 


(erg s"') 


(kms"') 


(kms-i) 


(IO^Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


3C120 


lwp04153 


44.46± 0.04 


2780 












lwp04500 


44.37± 0.03 


3837 












lwp05610 


44.35± 0.02 


2568 












lwp09048 


44.41± 0.03 


2638 












lwp09461 


44.20± 0.02 


3411 












lwp09850 


44.22± 0.02 


3399 












Iwp 10407 


44.29± 0.04 


4083 












Iwp 11524 


44.28± 0.03 


2828 












Iwp 11 946 


44.27± 0.03 


2786 












lwpl2536 


44.36± 0.03 


2412 












lwi01317 


43.86± 0.08 


2874 












lwr02983 


43.73± 0.02 


2759 












lwr06849 


44.12± 0.03 


3997 












lwr09102 


44.34± 0.03 


2875 












lwr09778 


43.93± 0.02 


2735 












lwrl3786 


44.30± 0.02 


2758 












lwrl5618 


44.25± 0.03 


3040 












lwrl6609 


44.17± 0.04 


3053 












lwrl6874 


44.39± 0.06 


4431 














44.23± 0.01 


3074 


44.09±0.09 


2327±50 


1166±50 




3C 390.3 


y33y0204t 


42.63± 0.17 


7884 


43.64 d=0. 14 


12694±13 


3105±81 


287 ± 64 


Aknl20 


y29e0305t 


44.48± 0.08 


4377 


43.93±0.04 


6143±42 


1921±60 


150 ± 19 


FairaU9 


yOya0104t 


44.30± 0.08 


3769 


43.94±0.10 


5999±66 


3787±196 


255 ± 56 


IC4329A 








42.89±0.15 


5964±134 


2476±226 


9.90t\lll 


Mrk79 


lwiO1320 


43.81± 0.02 


5057 










lwi06141 


43.60± 0.02 


4179 














43.71± 0.02 


4597 


43.65±0.03 


4858±38 


1882±121 


52.4 ±14.4 


MrkllO 


lwpl2760 


43.64± 0.08 


2216 












lwpl2761 


43.82± 0.07 


2504 














43.73± 0.05 


2355 


43.66±0.04 


1600±13 


1002±53 


25.1 ±6.1 


Mrk279 


lwp02522 


44.03± 0.03 


4812 












IwplOlie 


43.77± 0.02 


3314 












Iwp 15450 


43.03± 0.04 


5330 












lwpl5687 


43.75± 0.02 


3275 












lwpl9173 


43.98± 0.04 


3555 












lwpl9220 


44.02± 0.03 


3987 












lwpl9598 


44.28± 0.05 


3699 












lwpl9937 


44.20d= 0.04 


4204 












lwp2()271 


44.26d= 0.05 


4305 












lwp20725 


44.18±0.03 


5057 












lwr03073 


43.03± 0.10 


3710 












lwrl0816 


43.98± 0.04 


8201 












lwrll623 


43.94± 0.02 


6934 












lwrl5803 


43.96± 0.03 


4121 
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Table 1 — Continued 



Name 


DatalD 


logiaooo 


FWHM(Mgn) 


logLsioo 


FWHM(H/3, mean) 


o-/,„e(H/3, rms) 


logMBH(RM) 






(erg s"') 


(kms-l) 


(ergs-l) 


(kms-l) 


(km s-i) 


(10«Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 






43.89± 0.01 


4441 


43.66±0.08 


5354±32 


1420±96 


34.9 ±9.2 


Mrk335 


y29e0205t 


44.12± 0.04 


1977 












y29e0206t 


44.13± 0.05 


1783 














44.13± 0.03 


1878 


43.78±0.02 


1735±2 


933±64 


14.2 ±3.7 


Mrk509 


y0ya0305t 


44.55± 0.05 


3357 


44.16±0.10 


3015±2 


1276±28 


143 ± 12 


Mrk590 








43.46±0.05 


2906±89 


i '^'^ i 1 Ar\ 

1231 ±49 


47.5 ±7.4 


MrkolV 


Iwrl 1936 


44.07± 0.03 


3597 












lwrl3704 


44.03± 0.04 


4565 














44.05± 0.03 


4053 


/I O /'A 1 f\ i A 

43.64±0.14 


4899±37 


1680±67 


49.4 ± 7.7 


NuC3227 


o5kp01010 


41.76± 0.04 


3688 


42.48±0.04 


5103±160 


1925±124 


/(no 1 O 1 A 

42.2 ±21.4 


NGC3516 


y31iO105t 


43.05± 0.02 


3919 












y31iO206t 


43.23± 0.02 


3775 












y31r0306t 


43.19± 0.02 


3796 












y31r0406t 


42.73± 0.02 


3499 












y31r()5()6t 


43.08± 0.02 


3527 














43.06± 0.01 


3699 


42.62±0.28 


5840±1976 


1837±115 


42.7 ±14.6 


NGC 3783 


o57b01010 


43.39± 0.06 


2524 


43.02it0.06 


3770±68 


1753±141 


29.8 ±5.4 


NGC 4051 


IwplllOU 


41.76± 0.02 


1574 












lwpl2092 


41.70± 0.02 


1387 












Iwp 12092 


41.68± 0.02 


1503 












Iwp 19265 


41.58± 0.02 


1305 












lwp20497 


A 1 CI 1 1^ Ci'^ 

41.53± 0.02 


946 












lwp23153 


A 1 '~lf\ 1 C\ 

41.70± 0.02 


918 












lwp24347 


41.81± 0.02 


752 












lwp27297 


41.77± 0.02 


1349 












lwp2729» 


41.77± 0.02 


1667 












Iwr01728 


41.76± 0.02 


2582 














41.68± 0.01 


1322 


41.88±0.08 


654±2 


916±64 




NGC4151 


O42303070 


43.04± 0.03 


4905 










O59701040 


42.20± 0.02 


3020 














42.62± 0.02 


3849 


41.92±0.23 


6371±150 


1914±42 


45 7+^ '' 


NGC 4593 


lwp02731 


42.82± 0.02 


4123 










Iwp05348 


42.82± 0.02 


3022 












Iwp05371 


42.81 ±0.02 


3288 












Iwp05394 


42.83± 0.02 


2778 












lwp05411 


42.80± 0.02 


3444 












lwp05430 


42.74± 0.02 


3107 












lwp06266 


42.64± 0.02 


3237 












lwp06300 


42.72± 0.02 


3033 












Iwp 12278 


43.01 ± 0.02 


2822 












lwpl2279 


43.03± 0.02 


5194 












lwr07884 


42.83± 0.02 


2723 












lwr09818 


42.77± 0.02 


1929 












lwrl0539 


42.86± 0.02 


2987 
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Table 1 — Continued 



Name 


uataiiJ 


iOg 1.3000 


r W rlJVU^JVlg ii ) 


iOgI.5100 


rvvrijvi(_rip, mean) 


cr/me(np, rms) 


iOgMBH(KM) 






(erg ) 


A-TVl o"!^ 
^^KQl S ) 


(erg s~^) 


yKUl S } 


yKni S } 


(A Ci^M^ \ 
t^lU Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




Iwr 10622 


42.8 1± 0.02 


3117 












lwrl6177 


42.87± 0.02 


3321 














42.82± 0.01 


3140 


42.85±0.04 


5143±16 


1561+55 


9.8 + 2.1 


NGC5548 


y0ya0205t 


43.01± 0.03 


4756 


43.31±0.02 


6107±23 


2063+32 


65.4+?-6 

-2.5 


NGC7469 


y3b6010bt 


43.67± 0.03 


3061 


43.30±0.05 


1722±30 


1456+207 


12.2+1.4 


PG 0026+129 


y2jk0108t 


45.18± 0.04 


1104 


44.95±0.08 


2544±56 


1774+285 


393 + 96 


PG 0052+251 








44.78±0.12 


5008±73 


1783+86 


369 + 76 


PG 0804+761 


lwrl3645 


45.14± 0.02 


5175 












lwrl6666 


44.93± 0.12 


2533 














45.04± 0.06 


3621 


44.88±0.09 


3053±38 


1971+105 


693 + 83 


PG 0844+349 


yOp80105t 


44.53± 0.04 


3045 


44.19±0.07 


2694±58 


1448+79 


92.4 + 38.1 


PG 0953+414 








45.15±0.07 


3071 ±27 


1306+144 


276 + 59 


PG 121 1+143 


yOiz0403t 


44.81±0.03 


1610 












y0iz0404t 


44.82± 0.03 


1642 














44.82± 0.02 


1626 


44.70±0.08 


2012±37 


1080+102 


146 + 44 


PG 1226+023 


y0g4020et 


46.08± 0.06 


3420 












y0g4020ft 


45.98± 0.06 


3105 












y0g4020ht 


46.20± 0.09 


2900 












y0g4020jt 


46.24± 0.08 


2839 












y0g40201t 


46.24± 0.07 


3032 












y0g4020nt 


46.24± 0.07 


3022 












yOnb0104t 


45.94± 0.06 


3422 














46.13± 0.03 


3098 


45.93±0.06 


3509±36 


1777+150 


886 + 187 


PG 1229+204 


lwrl3136 


44.41± 0.07 


3054 












lwrl6071 


44.47± 0.03 


4940 














44.44± 0.04 


3884 


43.65±0.06 


3828±54 


1385+111 


73.2 + 35.2 


PG 1307+085 








44.82±0.05 


5059±133 


1820+122 


440 + 123 


PG 1411+442 


O65617010 


44.81± 0.03 


2452 


44 52+0 05 


2801 +43 


1 (507+ 168 


443 + 146 


PG 1426+015 


lwp05440 


45.14± 0.03 


6957 












lwp05446 


45.04± 0.06 


4575 












lwrl6020 


44.92± 0.07 


6776 














45.04± 0.03 


5997 


44.60±0.09 


7113+160 


3442+308 


1298 + 385 


PG 1613+658 


Iwp 19372 


45.09± 0.03 


7518 












lwpl9380 


45.10± 0.03 


7996 














45.09± 0.02 


7753 


44.73±0.10 


9074+103 


2547+342 


279 + 129 


PG 1617+175 


lwp07592 


44.81± 0.06 


5951 












lwp25629 


44.64± 0.05 


4375 














44.73± 0.04 


5102 


44.36±0.10 


6641+190 


2626+211 


594+138 


PG 1700+518 








45.56±0.03 


2252+85 


1700+123 


701+182 


PG 2130+099 


lwp02520 


44.60± 0.05 


3042 










lwp03568 


44.66± 0.08 


2208 












lwp07205 


44.53± 0.03 


2039 












lwr01774 


44.47± 0.06 


1820 












lwr04610 


44.52± 0.03 


1290 
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Table 1 — Continued 



Name 


DatalD 


logi-3000 


FWHM(Mgll) 


log L5 100 


FWHM(H/3, mean) 


criineW, rms) 


logMBH(RM) 






(erg s"') 


(km s-i) 


(erg s-i) 


(km s"') 


(km s"') 


(10«Mo) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




lwi04628 


44.58± 0.06 


2479 












lwrl5802 


44.51 ±0.03 


2885 














44.55± 0.02 


2174 


44.40±0.02 


2853±39 


1624±86 


38 ±15 



Note. — Column (1) object name; Column (2) identification name of the spectrum in lUE (prefix "Iw") or HST archives; Column (3) monochro- 
matic continuum luminosity at 3()()() A. For objects with multiple spectra, mean monochromatic continuum luminosity at 3000 A is listed in bold at 
the bottom of each object. Column (4) FWHM of broad Mg II. For objects with multiple spectra, mean FWHM is listed in bold at the bottom of each 
object. Column (5) monochromatic continuum limiinosity at 5100 A taken from Bentz et al. (2009). For objects having multiple measurements, 
here hsted is mean value. Column (6) H/3 FWHM of mean spectra taken from Collin et al. (2006) and recent update Usted in Section 3.2. For 
objects having multiple measurements, here listed is mean value. Colunm (7) H/3 <t;,„(, of rms spectra taken from Peterson et al. (2004) and recent 
update listed in Section 3.2. For objects having multiple measurements, here Usted is mean value. Column (8) BH mass from RM taken from 
Peterson et al. (2004) and recent update listed in Section 3.2. 
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Table 2. Continuum and Emission-line Parameters of the SDSS Sample 



SDSS Name 


z 




FWHM(H/3*) 


logF(H/3'') 


logF(H/3") 




FWHM(Mgn'') 


logFCMgn'') 


logF(Mgn») 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


JOOOOll.96+000225.3 


0.478 


44.69 


3037 


-13.89 


-15.68 


44.99 


3284 


-13.84 


-14.89 


JOOOl 10.97-105247.5 


0.529 


44.98 


6807 


-13.73 


-15.53 


45.18 


5797 


-13.85 


-15.15 


J001725.36+141132.6 


0.514 


45.23 


5676 


-13.49 


-15.52 


45.49 


4432 


-13.64 


-14.71 


J002019.22-1 10609.2 


0.492 


44.85 


2832 


-13.74 




45.00 


2677 


-13.91 


-16.16 


J005121.25+00452I.5 


0.727 


45.04 


2572 


-14.09 


-15.07 


45.14 


1606 


-14.42 


-15.41 


J005441. 19+0001 10.7 


0.646 


45.08 


2220 


-14.24 




45.13 


2172 


-14.20 


-15.29 


J010448.57-091013.0 


0.469 


44.77 


4610 


-13.81 


-15.33 


44.88 


2627 


-14.11 


-15.86 


J010644.16-103410.6 


0.468 


44.72 


3873 


-13.83 


-15.35 


44.85 


3074 


-13.80 


-1543 


J011132.34+133519.0 


0.576 


45.13 


8060 


-13.66 


-15.72 


45.39 


5495 


-13.73 


-14.87 


J012016.73-092028.8 


0.495 


44.71 


3284 


-13.72 


-16.02 


45.05 


3312 


-13.58 


-16.07 



Note. — Column (1) object name; Column (2) redshift derived from the peak of [OIII] A5007. Column (3) luminosity of the power-law continuum at 5100 A. Column (4) FWHM of broad H/5. 
Column (5) flux of the broad component of H/3. Column (6) flux of the narrow component of H/3. Column (7) luminosity of the power-law continuum at 3000 A. Column (8) FWHM of broad 
Mgll. Column (9) flux of the broad component of Mg II. Column (10) flux of the narrow component of Mg n. Luminosities, fluxes, and FWHM are in units of erg s"\ erg s""^ cm""^, and km s"^, 
respectively.(This table is available in its entirety in a machine-readable form in the onhne journal. A portion is shown here for guidance regarding its form and content.) 
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Table 3. Regression Results for 



log 



FWHM(Mgll) 
1000 km s-i 



= klog 



FWHM(H/?) 
1000 km s-i 



+ c 



Method 


k 








c 




OLS 


0.73 ±0 


02 





09 


±0 


01 


WLS 


0.73 ±0 


02 





09 


±0 


01 


FlTexy 


0.77 ±0 


01 





06 


±0 


02 


FITexy_T02 


0.81 ±0 


03 





04 


±0 


02 


Gaussfit 


0.81 ±0 


03 





04 


±0 


02 


BCES (bisector) 


0.78 ±0 


03 





06 


±0 


02 


BCES (orthogonal) 


0.81 ±0 


02 





05 


±0 


01 


LINMIX ERR 


0.79 ±0 


03 





05 


±0 


02 



Table 4. Regression Results for log 



Mbh(RM) 
10* Mq 

Mbh Comparisons 



a + piog 



^3000 



lO'*^ erg s- 



+ 7 log 



[1S¥?] and 



Scheme 


a 


P 


7 


AMbh(RM) 


AM 


bh(H/3) 


AMbh(H/3) 


AMbh(H/3) 














(Vestergaard+06) 


(Collin+06) 


(Ours) 


(1) 


(2) 


(3) 


(4) 


(5) 






(6) 


(7) 




(8) 










Mean 


1 cr 


Mean 


1 a 


Mean 


1 a 


Mean 1 a 


Scheme 1^ 


1.15 ±0.27 


0.46 ±0.08 


1.48 ±0.49 


0.01 


0.39 


0.15 


0.22 


0.10 


0.14 


0.07 0.12 


Scheme 2 


0.88 ±0.08 


0.48 ±0.08 


2 


0.01 


0.42 


0.11 


0.19 


0.06 


0.18 


0.03 0.18 


Scheme 3 


1.03 ±0.08 


0.48 ±0.08 


1.70 


0.01 


0.40 


0.13 


0.20 


0.07 


0.15 


0.05 0.14 


Scheme 4 


1.13 ±0.27 


0.5 


1.51 ±0.49 


0.01 


0.40 


0.10 


0.21 


0.05 


0.14 


0.03 0.12 



Note. — Fits for the 29 objects with Mgll data measured in the paper. AMbh = log Mbh -log Mbh (Mg II) are the differences of the masses 
obtained from different methods with masses estimated by our Mg II formalism for each scheme. Column (5) The mean and standard deviation of 
the deviations between BH masses obtained from RM and our single-epoch Mgll estimators. The mean and standard deviation of the deviations 
between masses estimated from our single-epoch Mg II estimator and masses derived using (Column 6) the H/3 formalism of Vestergaard & Peterson 
(2006), (Column 7) the H/3 formalism of Collin et al. (2006), and (Column 8) the new H/3 formahsm obtained in this work (Equation (11)) for the 
SDSS sample. 



"Scheme 1 is fitted by using the code MLINMIX_ERR of Kelly (2007). 



